Numerical cubature from Archimedes' hat-box theorem 
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Archimedes' hat-box theorem states that uniform measure on a sphere projects to uniform measure on an 
interval. This fact can be used to derive Simpson's rule. We present various constructions of, and lower bounds 
for, numerical cubature formulas using moment maps as a generalization of Archimedes' theorem. We realize 
some well-known cubature formulas on simplices as projections of spherical designs. We combine cubature 
formulas on simplices and tori to make new formulas on spheres. In particular S" admits a 7-cubature formula 
(sometimes a 7-design) with 0(n 4 ) points. We establish a local lower bound on the density of a PI cubature 
formula on a simplex using the moment map. 

Along the way we establish other quadrature and cubature results of independent interest. For each t, we 
construct a lattice trigonometric (2f + l)-cubature formula in n dimensions with 0(n') points. We derive a 
variant of the Moller lower bound using vector bundles. And we show that Gaussian quadrature is very sharply 
locally optimal among positive quadrature formulas. 



1. INTRODUCTION 

Let fi be a measure on M." with finite moments. A cubature 
formula of degree t for jj. is a set of points F — {p a } C K" and 
a weight function p a h-> w a £ K such that 



/ 



P{x)dfL = P(F) 



dcf 



Y, WaP(Pa) 



for polynomials P of degree at most t. (If n = 1, then F is also 
called a quadrature formula.) The formula F is equal-weight 
if all w a are equal; positive if all w a are positive; and negative 
if at least one w a is negative. Let X be the support of fx. The 
formula F is interior if every point p a is in the interior of X; 
it is boundary if every p a is in X and some p a is in dX; and 
otherwise it is exterior. We will mainly consider positive, in- 
terior (PI) and positive, boundary (PB) cubature formulas, and 
we will also assume that jj, is normalized so that total measure 
is 1. PI formulas are the most useful in numerical analysis 
I2H Ch. 1]. This application also motivates the main question 
of cubature formulas, which is to determine how many points 
are needed for a given formula and a given degree t . Equal- 
weight formulas that are either interior or boundary (EI or EB) 
are important for other applications, in which context they are 
also called t-designs. 

Our starting point is a connection between quadrature on 
the interval [—1,1] and cubature on the unit sphere S 2 , both 
with uniform measure. By Archimedes' hat-box theorem Q], 
the orthogonal projection % from S 2 to the z coordinate pre- 
serves normalized uniform measure. In plainer terms, for any 
interval I C [a,b] or other measurable set, the area of 7t (I) 
is proportional to the length of /; see Figure^ (It is called the 
hat-box theorem because the surface area of a hemispherical 
hat equals the area of the side of a cylindrical box containing 



it.) Therefore if F is a f-cubature formula on S 2 , its projection 
7t(F) is a f-cubature formula on [—1, 1]. 




Figure 1: Archimedes' hat-box theorem. 



The 2-sphere S 2 has 5 especially nice cubature formulas 
given by the vertices of the Platonic solids. Their cuba- 
ture properties follow purely from a symmetry argument of 
Sobolev 12511 - Suppose that G is the group of common sym- 
metries of a putative cubature formula F and its measure fl . If 
P(x) is a polynomial and Pg(x) is the average of its G-orbit, 
then 



P G (x)dn = / P{x)dfi P G {F) = P(F). 
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Therefore it suffices to check F for G-invariant polynomials. 
In particular, if every G-invariant polynomial of degree < t is 
constant, then any G-orbit is a f-design. 

By Sobolev's theorem, the vertices of a regular octahe- 
dron form a 3-design on S 2 . If we project this formula using 
Archimedes' theorem, the result is Simpson's rule. Another 
projection of the same 6 points yields 2-point Gauss-Legendre 
quadrature. Figure [2] shows both projections. The 8 vertices 
of a cube are also a 3-design. One projection is again 2-point 
Gauss-Legendre quadrature; another is Simpson's | rule. Fi- 
nally the 12 vertices of a regular icosahedron form a 5-design 



(a) Simpson's rule 



(b) 2-point Gauss-Legendre rule 



Figure 2: Two projections of the octahedron rule. 



by symmetry. One projection of these 12 points is 4-point 
Gauss-Lobatto quadrature. 



The rest of this article applies toric moment maps, which 
generalize Archimedes' theorem to higher dimensions, to the 
cubature problem. Section |2] shows that several well-known 
quadrature formulas on the interval and cubature formulas 
on simplices are projections of higher-dimensional, symmet- 
ric formulas. Section [4] combines formulas on tori with for- 
mulas on simplices and moment maps to make formulas on 
spheres and projective spaces. In particular, it constructs a PI 
7-cubature formula on the sphere S" with 0(n 4 ) points. Fi- 
nally Section |6] uses moment maps to establish a local lower 
bound for the density of points in any PI cubature formula on 
a simplex. A similar lower bound holds for an arbitrary simple 
convex polytope. 



Along the way we establish some other quadrature and cu- 
bature results that are not derived from moment maps but are 
of independent interest. Section[3]establishes new lattice cu- 
bature formulas on tori that are similar to cubature formulas 
based on error-correcting codes I17ll . In particular it con- 
structs, for each t, a trigonometric (2? + l)-cubature formula 
on [Q,2n) n of lattice type with 0(n') points. This improves 
a construction of Cools, Novak, and Ritter with 0(n 2 ') points 
and negative weights |5], and agrees up to a t -dependent con- 
stant factor with the Stroud-My sovskikh lower bound 1 191 12911 . 
Section [5] presents a refinement of this well-known lower 
bound in odd degree. It is similar to the Moller bound 1 18], 
but applies to some new cases. Section[6]also establishes that 
Gaussian quadrature is very sharply locally optimal among 
all positive quadrature formulas (Theorem l6.3l >. This bound 
might be previously known since Gaussian quadrature has 
been widely studied, but the author could not find it in the 
literature. 



2. PROJECTION CONSTRUCTIONS 

The immediate higher-dimensional generalization of 
Archimedes' theorem replaces the sphere S 2 by the complex 
manifold CP". This manifold has a natural metric and a nat- 
ural real algebraic structure. Concretely, assume that the pro- 
jective coordinates (zo '■ Zi : ••• : z n ) of CP" are normalized so 
that 

ko| 2 + ki| 2 + ... + k„| 2 = i. 

Then the coordinates ZkZj together embed CP" into C'" +1 '" 

as a real algebraic variety (with C^ +1 ' 2 interpreted as a 
2(n + l) 2 -dimensional real vector space) and a Riemannian 
manifold. This embedding is familiar in quantum mechan- 
ics as the density matrix (or density operator) formalism l20l 
§2.4]. The induced metric is called the Fubini-Study metric. 
Since the metric yields a measure on CP", and since it is a real 
algebraic variety, we can consider cubature formulas on it. 
There is a projection % : CP" — * A„ to the n-simplex given 

by 

n(zQ :zi : ••• :z n ) = (\zo\ 2 ,\zi\ 2 ,...,\z n \ 2 ), 

using normalized coordinates for CP" and barycentric coor- 
dinates for A„. It is linear and it preserves normalized mea- 
sure. In more abstract terms, % has these properties because 
CP" is a projective toric variety and % is its moment map. 
Archimedes' theorem is a description of the moment map of 
CP 1 = S 2 . Thus, if F is an interior f -cubature formula on CP", 
then 7C(F) is a f-cubature formula on A n . 

Ivanovic, Wootters, and Fields 1 12, 32] defined one interest- 
ing family of 2-designs on CP 9-1 for q = p k a prime power. If 
p is odd, then the 2-design is the orbit of a standard basis vec- 
tor ek in the group generated by cyclic permutation and linear 
operators of the form 

L(e k ) = co Tl " {ak2+bk+c) e kl 
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Figure 3: A 2-dimensional generalization of Simpson's rule. 



where co is a pih root of unity and Tr p is the ¥ p trace function 
on Fo. The construction is more complicated when p = 2. In 
either case, the standard basis projects to the vertices of A^_i 
and the other q 2 vectors project to the center. The result is a 
standard degree 2 generalization of Simpson's rule for A 9 _i, 
shown in Figure[3]when q = 3. 

Other interesting designs and cubature formulas on CP" - 1 
come from designs and formulas on S 2 "~ 1 . The generalized 
Hopf fibration 

h-.s 2 "- 1 ^cp"- 1 

is a quadratic, volume-preserving map from S 2 " -1 to CP" . 
Namely, if we place S 2 "~ l in C", h takes each point to the 
complex line containing it. The map h projects a It- or It + 1- 
cubature formula on S 2 " -1 to a t -cubature formula on CP" -1 . 

One interesting example is the 240 roots of the Eg root sys- 
tem, which are a 7-design as well as the solution to the sphere 
kissing problem in R 8 A §14.2]. The root system has two 
natural positions in C 4 . In the first position, it is generated 
from the two points 

(1,1,1,0) (l-©,0,0,0) 

by freely permuting the first 3 coordinates, applying the map 
(a,b,c,d) h- > (d,a,—b,c), and multiplying any one coordinate 
by co, a cube root of unity. In the second position, it is gener- 
ated from the three points 

(1,1,1,1) (2,0,0,0) (l+/,l+/,0,0) 

by freely permuting the four coordinates and multiplying any 
two coordinates by i. These two positions respectively exhibit 
the Eisenstein and Gaussian lattice structures of the Eg lattice. 
The Hopf fibration sends the Eisenstein position of the root 
system to a 40-point 3-design in CP 3 and the Gaussian po- 
sition to a 60-point 3-design. Then the moment map projects 
these two 3-designs to 3-cubature formulas for the tetrahedron 
A3 that appear in Abramowitz and Stegun (ji p. 895]. They 
have 8 and 1 1 points, respectively, and are shown in Figure^ 
The composition n oh of the moment map and the Hopf 
fibration is a torus fibration T 2 : R 2 " — > A„_ 1 that does not fully 
depend on the complex structure R 2 " = C", but only on the 
decomposition of R 2 " into n orthogonal planes. Explicitly, 
the map is: 

Ti(jci,. . . ■X2n) — {x\ + x\,x\ . . . ,x\ n _ Y +x\ n ). 



This projection is analogous to a map Tj : S" — > A„_i de- 
fined by Xu |H : 

Tl[xi, . . . ,X n ) = (x 2 , . . . ,X 2 ). 

The Xu map does not preserve uniform measure. Rather, 
it takes uniform measure on the sphere to the measure with 
weight function 

2 n n n/2 

in barycentric coordinates. 

In the case of the Eg root system, one interesting set of or- 
thogonal planes is the 4 eigenplanes of the abelian subgroup 
of Aut(£g) of the form C$ x C5. Eric Rains 12311 has computed 
the corresponding 3-cubature formula on A3 using Magma 
ll33ll . In barycentric coordinates on A3, its points and weights 
are the orbits of the two weighted points 

£ 1= 1(0,0,5 -V5,5 + V5) Wl = ^ 

P2 = 1(2,2,3 + v/5,3 -y/5) w 2 = l 

under the action of the coordinate permutations (34) and 
(13)(24). In particular, it has 8 points. In conclusion, at least 
three interesting 3-cubature formulas for A3 arise as projec- 
tions of Eg root system. The root system model explains the 
simple rational values of the weights. 

The Eg lattice is one of four widely studied and highly 
symmetric lattices in low dimensions; the other three are the 
Coxeter-Todd lattice K12 in R 12 , the Barnes-Wall lattice Ai6 
in R 16 , and the Leech lattice A 24 in R 24 |4, Ch. 4]. In each 
case, cases, the set of short vectors has transitive symmetry, 
and in each case, Sobolev's theorem establishes its degree as 
a spherical design. 

The 756 short vectors of K\2 form a 7-design on S x 1 . In one 
of its several presentations as an Eisenstein lattice in C 6 (the 
"3-base" presentation |4, §7.8]), the short vectors are gener- 
ated from the two points 

(1,1,1,1,1,1) (1 -co,co- 1,0,0,0,0) 

by freely permuting coordinates, multiplying the coordinates 
by powers of co whose exponents sum to 0, and negating all 
coordinates. The projection T2 sends these points to a 16-point 
3-cubature formula on A5 generated from the points 

Pi = i(l,l, 0,0,0,0) W1 = J2 

p 2 = ^(l,l,l,l,l,l) w 2 = g 

by freely permuting coordinates. This formula was found by 
Stroud ll27ll28ll. 

The 4320 short vectors of Aig form a 7-design on S 15 . In 
its simplest position (which exhibits its Gaussian lattice struc- 
ture), the short vectors are generated from the two vectors 

(1,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0) 
(2,2,0,0,0,0,0,0,0,0,0,0,0,0,0,0) 



(a) 8 points 



(b) 11 points 



Figure 4: 3-cubature formulas on A3 from the Eg root system. 



by permuting coordinates under the group GL(4,2) k (Z/2) 4 
of affine automorphisms of (Z/2) 4 , together with sign 
changes that keep the coordinate sums divisible by 4. The 
projection T2 sends these points to a 51-point 3-cubature for- 
mula on A7 generated from the points 



Pi = 


(1,0,0,0,0,0,0,0) 


W\ 


P2 = 


1(1,1,0,0,0,0,0,0) 


w 2 


P3 = 


1(1,1,1,1,0,0,0,0) 




PA = 


1(1,1,1,1,1,1,1,1) 


Vi'4 



1 



1080 
1 

270 
4 

135 
64 

135' 



under the action of the affine group GL(3,2) x (Z/2) 3 . This 
is not an optimal PI 3-cubature formula, because the orbit of 
P2 can be eliminated, leaving only 23 points. But it does have 
a novel property: Instead of full symmetrization, the orbit of 
P3 is in the pattern of the (8,4,3) Steiner system. But this is 
as good as full symmetrization for 3-cubature, because any 
monomial of degree 3 involves at most 3 coordinates. The 
structure of this Barnes-Wall projection led the author to re- 
late cubature to combinatorial f -designs and orthogonal arrays 
0. 



The above position of Ajg is compatible with its Gaussian 
lattice structure. Eric Rains found another interesting position 
which is compatible with an Eisenstein lattice structure. The 
corresponding 3-cubature formula on A7 has 50 points. They 



are generated from 

1 

720 
1 

90 
1 

80 
1 

60 
1 

40 
1 

30 
1 

30 

by the coordinate permutations (12), (13) (24) (57) (68), and 
(15)(26)(37)(48). 

The 196560 short vectors of the Leech lattice form a 11- 
design on S 2i . The lattice has a space Eisenstein lattice struc- 
ture which Conway and Sloane call the complex Leech lat- 
tice §7.8]. The complex basis that they give leads to a 
5-cubature formula on An generated by the points 



Pi = 


(1,0,0,0,0,0,0,0) 


W\ 


P2 = 


1(1,1,1,1,0,0,0,0) 


W>2 


P3 = 


1(1,1,0,0,1,0,0,0) 


Vl'3 


PA = 


^(4,0,4,0,1,1,3,3) 


H'4 


P5 = 


^(4,0,4,0,1,1,1,1) 


Vl'5 


P6 = 


1(3,1,3,1,1,1,1,1) 


U'6 


Pi = 


^(3,1,1,1,1,1,3,1) 


W-j 







Pi = 


5- 






P2 = 


^< 




1 


P3 = 


78 




1 


PA = 


18 



1 



(4,4,1,1,1,1,1,1,1,1,1,1) 



W\ = 

W2 = 
W3 = 
W4 = 



10920 
9 

3640 
27 

1820 
27 

3640 



by the action of the Mathieu group M\2- In other words, the 
coordinates of p2 are permuted in the pattern of the (12,6,5) 
Steiner system and the points of the other coordinate are per- 
muted freely. The total is 276 points. Another interesting ba- 
sis of plane consists of the mutual eigenplanes of the (Z/5) 3 
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subgroup of the isometry group of the Leech lattice. Eric 
Rains has computed that the corresponding 5-cubature for- 
mula on An has 498 points, consisting of 22 orbits of the 
surviving coordinate permutations. However, since none of 
the Barnes-Wall formulas on A7 are optimal, it is not clear 
that the smaller of these formulas is either. 



3. TORUS CONSTRUCTIONS 

The constructions in the next section depend on an auxil- 
iary case that generally works out better than cubature on S", 
CP", and A„, namely cubature on algebraic tori. There is a 
developed theory for a special case of this problem known as 
trigonometric cubature |6, 7]. We will describe a more general 
class of problems, with one new result for the classic trigono- 
metric cubature problem (Theorem l3.3> . 

Consider a torus group T = (S 1 )" together with a faithful 
linear action on some real vector space V = Mr. Then we can 
identify T with any faithful orbit O to give it a real algebraic 
structure. Since T is a compact group, it also comes with 
Haar measure (i.e., uniform measure). Given both structures, 
we can then consider cubature formulas for T . If a cubature 
formula F is a f-design and forms a subgroup of T, then it is 
called a lattice formula, or an additive t-design. 

Proposition 3.1. The lattice cubature problem on T is equiv- 
alent to a lattice packing problem as follows: 

1. The real algebraic structure on T does not depend on the 

orbit O or the base point chosen on O. The ring of poly- 
nomials on T is the same as the character ring R(T). 

2. Every character % : T — ► C is homogeneous as a polyno- 

mial on T. Its degree defines a norm on T, the character 
group of T. The norm is generated by unit steps corre- 
sponding to the characters that appear in V (g> C. 

3. The characters that are constant on a subgroup F C T form 

a sublattice F C.T. This correspondence is a bijection 
between finite subgroups and sublattices such that \F\ = 
[f : F]. 

4. The subgroup F is a t-design if and only ifF has minimum 

distance d = t + 1. 

The proof of Proposition B.ll is lengthy but routine and can 
be left as an exercise for the reader. It is essentially estab- 
lished in the literature when T = T(SO(2n)) acts on R 2n by 
separate rotations in n orthogonal planes. This case is equiv- 
alent to the (cubic) trigonometric cubature problem, defined 
as cubature formulas on the n-cube [0,27zQ" which are exact 
for trigonometric polynomials of degree t |7]. All of the argu- 
ments generalize without change. 

When T = 7/(SO(2«)), f is naturally identified with 1 n , 
and its norm is the l\ or taxicab norm. Another torus of in- 
terest to us is T = r(PSU(« + 1)), the group of diagonal uni- 
tary matrices with determinant 1 modulo its center. It acts on 
> interpreted as the space of in + 1) x (n + 1) complex 



matrices, by conjugation. In this case T = A„, the root lattice 
of PS U(« + 1), and its norm is defined by taking the roots of 
A n as unit steps. 

Theorem 3.2. Given a real algebraic torus T of dimension n, 
let K C T ® M be the real convex hull of the unit steps in T . 
Let 8i(K) be the lattice packing density ofK, and let Vol A' be 
the volume ofK normalized by T. Let t > and let d = t + 1. 
Then the best additive t-design F on T has at least 

d»(VolK) d>'(VolK) 1 

points. 

Theorem l3.2l has been noted independently by several peo- 
ple for trigonometric cubature, but may originally be due to 
Frolov |T[:| . In outline, a lattice F C T with minimum dis- 
tance d produces a packing of the dilated body ^K. The pack- 
ing density 8l (K) then yields a lower bound on the index of 
F. On the other hand, if A is the center lattice of the best pack- 
ing of K, then when t is large, ^A can be approximated by a 

sublattice of T. This establishes the upper bound. 

Note also that the best A has rational coordinates relative to 
T (or they can be made rational if A is not unique), because 
K is a rational polytope. Thus there exist special distances d 
such that the best F has exactly 

d n (Vo\K) 
2"8 L (K) 

points. Also if some d achieves exactitude, then so does kd 
for every k > 1 . 

If T = r(SO(2n)) is the standard cubic n-torus, then K is 
the n-cross polytope C*. For example, Minkowski established 
that the lattice packing density of the regular octahedron C3 is 
j|. So there exists an additive 5-design on r(SO(6)) with 38 
points fill El . 

Since C\ is a square, its packing density is 1 . Noskov El 
found the best discrete approximation to this packing for every 
distance d to obtain lattice rules for T(SO(4)). If d — 2s, then 
the best approximation is exact and there is a (2s — l)-design 
with 2s 2 points. If d = 2s + 1, then the best approximation 
corresponds to the tiling of 1? by the discrete l\ ball of radius 
s, or the tiling of the plane R 2 by certain Aztec diamonds, 
as shown in Figure |5ja). The ball and the corresponding 2s- 
design have s 2 + (s + 1 ) 2 points. 

Noskov's designs have a counterpart for r(PSU(3)), where 

r(PSU(3)) = A2 is the triangular lattice. If we identify A2 
with the Eisenstein integers Z[fl)], then the highest-density lat- 
tice with minimum distance d is the ideal generated by 

d d+l. 

When d = 2s, the dual (d — 1) -design has 3s 2 points and 
exactly matches the tiling of the plane by regular hexagons. 
When d = 2s + 1, it has 3s 2 + 3s + 1 points and corresponds 
to a tiling the plane by the hexagonal polyhex of order s (an 
"Afghan hexagon"), as shown in Figure|3b). 
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Theorem 3.3. Let t > 0. The torus T(SO(2n)) has a (2t + 1)- 
design with 0{nf) points. More precisely it has a 2t -design 
with (2«y (1 +o(l)) points as n —> °° and a (2t + \ )-design 
with twice as many points. The torus r(PSU(« + 1)) has a 
t-design with n'{\ +o(l)) points as n — > oo. 



We define the lattice 



Remark. Theorem l3.3l can be compared with a prior result by 
Cools, Novak, and Ritter [5], who obtained NI f-cubature for- 
mulas with 0(n') points and negative weights. Another com- 
parison is with the lower bound due to Stroud and Mysovskikh 
1119112911 for trigonometric 2f-cubature: 



F > 



(2«)'(l-o(l)) 



tl 



The Moller bound applies to trigonometric (2f + l)-cubature 
in its interpretation as cubature on T(SO(2n)) because it is a 
centrally symmetric algebraic variety. It yields: 



|F|> 



2(2«)'(1- (1)) 



Section|5]establishes an analogous lower bound for f-cubature 
on J(PSU(n)) (CorollaryES: 



"'(1-^(1)) 
\t/2V.[t/2\! 



F > 



Thus for each t, Theorem l3.3l is asymptotically optimal to 
within a constant factor, even though the lower bounds do not 
require F to be positive or interior. 

Proof. By Proposition l3.ll our task is to find suitable lattices 

in Z" = r(SO(2n)) and An = r(PSU(n + 1)). Our task is 
fulfilled by Craig lattices |4, §8.6] in A n and skew analogues 
of Craig lattices in Z". We describe the A„ case first. 

We can model A„ as the set of points in Z" +1 with zero co- 
ordinate sum. Let p > n + 1 be prime, and index the standard 
basis {e a } of Z" +1 by some subset N C Z/p. Let p > t > 0, 
and define a linear map : Z" +1 -> (Z/p)' by 

0(e fl ) = (a, a 2 ,..., a'). 



A w (A n ) = ker DA n . 

Plainly the index of A^(A„) is at most p' = n*(l +o(l)). (If 
n is large and p 1=3 n, it is p', because any lower power of p 
would violate the Stroud-Mysovskikh bound.) 

We claim that the distance of A^'(A„) is t + 1. To show 

this, we will show that <p is injective on the simplex A^ C 
Z" q 1 of non-negative vectors with coordinate sum t. A vector 
xGA n with root-step length at most t can be expressed as the 

difference of two vectors in A^; therefore injectivity shows 
that none of these vectors lie in A^ (A„). 

We can interpret a vector.? € Ai'' as a multiset of S over the 
set {0,...,rc} with |5| =t: if 



m a e 



then m a is the multiplicity of a 6 S. In this interpretation, <j> (x) 
is the list of power sums 

for 1 < k < t. By standard inversion formulas 1 26], these 
power sums determine the elementary symmetric functions of 
the elements of S when p > t, which are the coefficients of the 
polynomial 



aes 

Thus <p (x) determines S as a multiset and the vector x, and it 

is injective on A^ . 

For Z" (with the l\ norm), let p > 2n be prime. Index the 
standard basis {e a } of Z" by some subset N C Z/p such that 
N is disjoint from —N. Define <j) : Z" — * {Z/p)' by 



(p(e a ) = (a,fl 3 ,fl 5 



,2<-l 



), 



7 



and define 

A W (Z") = ker</>. 

Then the index of A^(Z„) is again at most (and usually ex- 
actly) p l = n (1 + o(l)). Its distance property can be ex- 
plained by embedding Z" isometrically into A2„ using the map 

a : e a h-> e a -e- a . 

Then 

AW(Z") = a- 1 (A( 2 ')(A„)). 

Since A (2 ')(A„) has distance at least 2f + 1, so does AW(Z n ). 
We can boost the distance to 2t + 2 by passing to its even-sum 
sublattice. □ 




r 



r 



Figure 6: A Hamming-like "plus" tiling. 



Remark. When t = 1, the number in the proof of The- 
orem |^] need not be prime, and the lattices AW(Z n ) and 
AW(A„) produce lattice tilings of the ball of l\ -radius 1 in 
Z™ and the combinatorial simplex A^ in A„. For exam- 
ple, when n = 2, they are equivalent to familiar tilings of the 
plus pentomino (Figure |6} and the triangle trihex. The plus 
tiling resembles combinatorial tilings coming from Hamming 
codes |4, §3.2]. More generally, Craig lattices resemble low- 
distance BCH codes. This resemblance is what led the author 
to Theorem l3.3l 



4. FIB RATION CONSTRUCTIONS 

The projection construction in Section[2] while instructive, 
is backwards in a sense: It is harder to make f-cubature for- 
mulas for CP" -1 and S 2n ~ l than for A„_i for most values of 
n and t. In this section we will use the same projections to 
lift cubature formulas to spheres and projective spaces from 
simplices. The construction also requires the definition and 
constructions of cubature formulas on tori from Section[3] 

Theorem 4.1. Let a : X — > Y be one of the three projections 
h, it, or X% and let T be a generic fiber. Let s = 2t + 1 



when X — S 2 "~ l and s — t when X = CP" -1 . Given an in- 
terior (or boundary) t -cubature formula F for Y and an in- 
terior s-cubature formula Ft for T, there is a twisted prod- 
uct s-cubature formula Fx = Ft k Fy for X. It satisfies 
\Fx I = \Ft I |Fy I and it inherits positivity from its factors. In 
the boundary case, \Fx \ < \Ft\ \Fy\. 

Note that in the three cases, T is isomorphic to S , 
r(SO(2«)), and r(PSU(n)), respectively. 

Proof. Let ay be the discrete measure on Y corresponding to 
the cubature formula Fy, and let Ox — a* (<7y ) be the pull-back 
of Oy to X. In other words, for each point p of weight w in Fy, 
Ox has a term consisting of uniform measure on the torus fiber 
a~ l (p). Also let fix and fly be normalized uniform measure 
on X and Y. 
We claim that 

/ P(x)dHx = I P(x)do x 
Jx Jx 

for any polynomial of P of degree s; in other words jXx and 
Ox are s-cubature equivalent 1 17]. If we assume the natural 
group structure on T, then it acts on X in each of the three 
cases with Y as the set of orbits. Then 

/ P(x)do x = / P T (x)do x / P(x)dix x = / P T (x)dHx: 
Jx Jx Jx Jx 

where Pt is the average of P with respect to the action of T. 
The polynomial Pt then descends to a polynomial Py on Y of 
degree t , and 

/ P T (x)dox=P(F Y ) f P T {x)dii x = f P Y (y)dii Y 

because a preserves measure. 

The measure Ox evidently has a twisted product s-cubature 
formula Fx = Ft k Fy given by replacing each fiber by a copy 
of Ft- (A singular fiber corresponding to a boundary point of 
T can be replaced by a projection of Ft.) Since pix and Ox 
are s-cubature equivalent, Fx is a cubature formula for fix as 
well. □ 

Remark. The proof of Theorem l4.1l is analogous to Sobolev's 
theorem with the finite group G replaced by the torus T. In- 
deed the argument works for any compact group. 

The simplest case of Theorem l4.1l is the Hopf map h. In this 
case the theorem says that a f-cubature formula F on CP"~ 
lifts to a (2? + l)-cubature formula F' on S 2n ~ l with (2t + 
2)\F\ points. This relation was also observed by Konig 1 15]. 

Corollary 4.2. The n-sphere S" has a 7 -cubature formula with 
0(n 4 ) points for all n, more precisely 4n 4 (l +o(l)) points. 
The 3-sphere S 3 has a (2s + 1)- cubature formula with 

f(s+l)(s 2 + 3) sodd 
\(s+ l)(s 2 + s + 2) seven 

points. 
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Proof. The simplex A„ has a 3-cubature formula with 0(n) 
points 1 17] constructed using Hadamard designs. This can be 
combined with the 7-design on T(SO(2n)) with <9(n 3 ) points 
provided by Theorem 13. 3 1 for a total of 0(n 4 ) points. More 
precisely, the formula on A„ has points at the corners, each of 
which lifts to 0(n) points; a point in the center, which lifts 
to (9(« 3 ) points; and 2n + o(n) points on [n/2\ -dimensional 
faces, each of which lift to 2n 3 (l +o(l)) points. Only the last 
family of points is significant and it comprises 4n 4 (l +o(l)) 
points. 

Noskov's formulas from Section|3]include a (2s + l)-design 
on the square torus T(SO(4)) with 2(s + l) 2 points. When 
s is odd, this can be combined with the Gauss-Lobatto s- 
quadrature formula on the interval A 2 with points. Two 
of the fibers are circles and can be replaced by 2(s + 1) points 
instead of 2(s + l) 2 points. The total is then (s + i)(s 2 + 3) 
points. When s is even, it can be combined with the Gauss- 
Radau ^-quadrature formula with ^ points. In this case one 
fiber is a circle. □ 

Remark. The first part of Corollary 14.21 actually yields a 
7-design on S"~ l with 0(n 6 ) points whenever there is a 
Hadamard matrix of order n. In this case the weights of the 



3-cubature formulas on A„_i are 



2(n+lj(n+2) 



at the faces, and 



An 

(«+l)(«+2) 



n(n+\)(n+2) 



at the corners, 



at the center. Thus the 



weights are all commensurable up to a factor of 2n 2 (note that 
n is even) and the cubature formula can be interpreted as a 7- 
design with this multiplicity factor. Moreover, copies of the 
lattice formulas on the torus fibers can be shifted to make the 
design multiplicity-free. Better yet, the design need only have 
0(n 4 ) points if, for example n — A-l k . In this case the prime 
p used in the proof of Theorem 13. 31 can be replaced by the 
prime power l k+ . The number of points on each fiber then 
compensates for all but a bounded part of the factor of 2n 2 in 
the weights. 

The previous best construction of 7-designs on S is due 
to Sidelnikov l24ll and requires 0(2 ki > k+l ^ 2 ) points when n = 
2 k . 

A useful variant of Theorem l4. ll involves the moment map 
T2 : R 2 " — > R" defined by the same formula as x 2 ■ S 2 "~ l — > 
A„_i, namely: 



T 2 (xi,...,X 2 n) = {x{ + x\,x\ + 



14, ■ 



AnY 



This %2 takes uniform measure on the ball B„ to uniform mea- 
sure on the simplex 



{x\x k > 0,Y, x k < 1 



}• 



When n = 2, Noskov's formulas together with some ad hoc 
cubature formulas for the triangle yield some economical for- 
mulas for the 4-ball B4. For example, there is a PB 3-cubature 
formula on the triangle x,y > 0,x+y < 1 with points and 
weights generated from 



- I 2 2 ^ 

_ .161 + nVu 

P2 = { 1344 ' 0) 



25 

W] = — 
1 48 



W2 



16-2^14 
25 



by switching the coordinates and negating v 14. This formula 
lifts to 1 generic fiber in B4 which can be replaced with 32 
points and 4 singular fibers which are circles and can be re- 
placed with 8 points each. The result is a PI 7-cubature for- 
mula on B4 with 64 points. 

Wandzura and Xiao ll3lll found competitive PI s-cubature 
formulas for s up to 30; Figure shows one example. Most 
of these yield competitive PI (2s + 1) -cubature formulas on 
B4 and S 5 . The formulas could probably be improved further 
with a search on the triangle that favors nodes on the edges. 

The map T2 : K 2 " — > M" also takes Gaussian measure on M 2 " 
to exponential measure on W\_. For example, there is a PB ex- 
ponential 4-cubature formula on Ml 2 ^ with points and weights 
generated from 

pi w (1.50766353, 1.50766353) w x k, 0.354104443 

p 2 « (6.29508677, 1.76717584) w 2 w 0.00876905581 

p 3 w (0.285606152,0) w 3 ~ 0.556110610 

P4 w (3.27491992,0) w 4 ~ 0.0722468398 

by switching the coordinates. It lifts to 3 generic fibers with 
50 points each and 4 singular fibers with 10 points each. The 
result is a positive Gaussian 7-cubature formula on R 4 with 
190 points. 



5. AN ALGEBRAIC LOWER BOUND 

Let X be the Zariski closure of the support of a measure 
/X on W and let A be the ring of polynomial functions on 
X. (Recall that the Zariski closure of a set, or closure in the 
Zariski topology, is the smallest algebraic variety containing 
it.) In other words, A is the quotient of M.[x] by the ideal Ix 
of polynomials that vanish on X. The ring A has a degree 
filtration coming from the degree filtration of E.[x}. Stroud 
I28U29I1 established an important lower bound on an arbitrary 
2t -cubature formula F for jx (not necessarily positive or inte- 
rior): 



Theorem 5.1 (Stroud). 

then 



If F is a 2t -cubature formula for ju, 



IF I >dimA 



<t- 



Mysovskikh 1 19] observed that this applies to trigonometric 
cubature by taking X = T(SO(2n)). (And according to Moller 
lH8ll . the bound was noted independently in special cases by 
other authors, e.g., Radon.) 

Proof. Define a bilinear form 

b:A< t xA< t ^R 

by 

b(P,Q)= [ P(x)Q(x)d i i. 
Jx 
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The form b is positive-definite because the integrand of b(P, P) 
is non-negative; moreover if the integrand vanishes on X, then 
P = as an element of A. Therefore b is non-degenerate, and 
its rank is dimA< r . On the other hand, the integrand lies in 
A<2r, so a 2f-cubature formula F leads to the formula 



b{P,Q)=^kP(pk)Q{Pk)- 



This formula realizes I? as a sum of \F\ rank 1 forms. There- 
fore \F\ is at least the rank of b, as desired. □ 

An interesting scholium of the proof of Theorem l5.1l is that 
if F is a 2f-cubature formula, then its points suffice to interpo- 
late polynomials on X of degree t . 

It is curiously difficult to improve the Stroud bound for odd- 
degree cubature. However, the inference that lower bounds 
improve mainly in even degrees is not consistent with the 
Hopf fibrations 



Proof. The vector space M (which is naturally an A -module) 
has an A -valued Hermitian inner product a induced by the 
Hermitian inner product on V. More precisely, let V be the 
representation conjugate to V and let 

M = Inv G (y<8>B) 

be the corresponding conjugate of M. (Note that A and B are 
both self-conjugate by hypothesis.) Let 

e:V®V — >C 

be the linearization of the standard Hermitian inner product 
on V, and let 



m:B®B- 



B 



be the linearization of multiplication on B. Let a 1 be the com- 
position 



B<®V<&V( 



)B '^B®B- 



B. 



h:S 



2n+l 



CP" h : r(SO(2n + 2)) -> T(PSU(« + 1)). We can restrict the domain to 



On the one hand, these maps are quadratic and double the de- 
gree of cubature in passing from the target to the domain; in 
particular, they do not preserve odd and even. On the other 
hand, Sections|2]and|4]together show that cubature in the do- 
main and target are comparably difficult when n^> t. 

The Hopf fibration example suggests a generalization of 
Stroud's theorem involving group actions and degree dou- 
bling. 

Theorem 5.2. Let pi be measure on W, and let X be the 
Zariski closure of its support. Let Y C M. k be another affine 
real algebraic variety on which a compact group G acts. Let 
A and B be the rings of complex-valued polynomials on X and 
Y, and suppose that there is a ring isomorphism 

A Inv G (B) 

that doubles the filtration degree of A. Let V be a unitary 
representation of G and define the filtered vector space 

M = lm G (B®V). 

IfF is a t-cubature formula for p, then 

dimM< r 



F > 



dimV 



We will always take Y to be the coordinate ring of another 
algebraic variety Y which is a principal G-bundle overX, such 
that the bundle projection a : Y — > X is quadratic. The A- 
module M can then be understood as the space of polynomial 
sections of a vector bundle E over X with fiber V. The sec- 
tions in M<, then behave like polynomials elements of A, ex- 
cept that their degrees are half-integers. If Y = X and G is 
trivial, then E is the trivial line bundle and Theorem |^| re- 
duces to Theorem 15. II The hypotheses of Theorem l5.2l have 
been chosen so that the proof of Theorem 15.11 generalizes to 
the case when E is not trivial. 



M(g>M = Inv G (B <g> V ) ® Inv G (V ® B) . 

Since the restricted domain is G-invariant, we can then restrict 
the target to A. Let a be this restriction of a' . Although given 
as a linear map on M ®M, it can be reinterpreted as a Her- 
mitian inner product on M. In more geometric terms, if M 
comes from a bundle E overX with fiber V, then a(f,g) is the 
pointwise inner product of two sections / and g of E. 
Note that a is positive-definite in the sense that 

«(/,/)(*) >0 

for all x G X, and if a(fj) = 0, then / = G M. The rest 
of the proof follows that of Theorem 15. II Define a complex- 
valued Hermitian inner product b on M< t by 



b{f,g) = / a(f{x),g{x))dll. 

Then b is also positive-definite, because a is positive-definite 
and n is Zariski-dense in X. Thus, b has rank dimM< r . A 
cubature formula F realizes b as a sum of \F\ terms of rank at 
most dimV. □ 

We state three special cases of Theorem l5.2l as corollaries: 

Corollary 5.3. If \F\ is a (2t + 1)- cubature formula on CP", 
then 



n + t\ fn + t+1 



F < 



Proof. Let Y be S 2n+l and let G to S 1 C C acting by complex 
rotation on C" +1 D S 2n+1 . The bundle projection a : Y — > X is 
the Hopf map h. Let V = L\ be the tautological representation 
of S , so that E is the tautological line bundle on CP". The 
space M<2 r +i is explicitly realized as the space of homoge- 
neous polynomials in z and z of bidegree (t+ l,t). The result 
follows by noting that 



dimM. 



<2t+l 



n + t+l 
n 



and that dimLi = 1. 



□ 
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Corollary 5.4. If \F\ is a t-cubature formula on r(PSU(n 
1)), then 



\F\> 



\t/2\\\t/2\\ 



Proof. Let Y be the torus T(SO(2n + 2)) and let G = S l again 
act by complex rotation in C" . In this case M<2t+i is 
spanned by the space of monomials in z and z of bidegree 
(s + 1 , s) with s < t and with the relation 

ZkZk = 1 

for all k. Its dimension is the number of points in the 
Minkowski difference 

A ('+i)_ A W 

where A„ is the discrete simplex defined in the proof of The- 
orem l3.3l This is very similar to Theorem l5.1l for 2f-cubature, 
because 



dimA 



<t 



|A«-A«| 



There is no concise formula for either number, but there is a 
concise estimate for fixed t in the limit n — > °°. If £ is either 
the trivial bundle when t is even or the bundle L\ (restricted 
from CP") when t is odd, then 



dimM< 



n+ 1 



St/2\,\t/2\,n + l-t) \t/2V.[t/2\l 
as t — > oo, as desired. □ 



Corollary 5.5. // /I is a Zariski- dense measure on S" and F 
is a (2? + \ )-cubature formula for ju, then 



Remark. When F is a lattice formula, Corollary 15. 4l is equiv- 
alent to Minkowski's classic upper bound on the density F as 
a lattice packing of the discrete simplex A^. This and the 
fact that the Hopf fibration is quadratic led the author to The- 
oreml5.2l 



|F|>2 



n—l+t 



Proof. The idea is to let Y = Spin(« + 1) and G = Spin(n), 
where Spin(n) is the connected Lie group that double covers 
SO(n). Then 

X = Spin(« + 1 ) /Spin(n) = SO(n + 1 )/SO(n) = S" . 

Our choice for the representation V is the spinor representa- 
tion of Spin(«) when n is odd and a semispinor representation 
when n is even. The rest of the argument is a review of stan- 
dard but non-trivial representation theory 1 14l I30B : 



1. The irreducible (unitary) representations of the spin groups 

Spin (2/1 ) and Spin(2n + 1) are both indexed by vectors 

A = a i X\ + «2 A2 + . . . + a„A„ 

with each a^ a non-negative integer. The representation 
with highest weight A can be called Vq(X) with G = 
Spin (2/1 ) or G = Spin(2n + 1). 

2. The representation Vc(tX\) is realized as the space A t 

of homogeneous polynomials of degree t on S 2 "^ 1 or 

S 2n . The representation Vspin(2«+i)(^n) i s tne spinor 
representation of dimension 2". The representations 
Vs P in(2n)(A n -i) andy Spin(2 „ ) (A„) are the semispinor rep- 
resentations of dimension 2"~ 1 . 

3. For any choice of fiber V — Vspin(n) (A), the A-module M is 

the induced representation 

M = M(n,A)=Indg; +1 V Spin(n) (A). 

4. The characters of V (A ) for all A generate the natural real 

algebraic structure on Spin(2n) or Spin(2« + 1) (indeed 
on any compact Lie group). One degree filtration is de- 
fined by letting the characters of every V(X/ C ) have de- 
gree 2, except for the spinor or semispinor representa- 
tions, which have degree 1 . 

5. The structure of M(«,A) can be computed by branching 

formulas for the restriction of an irrep of Spin(« + 1) 
to Spin(n), together with induction-restriction duality. 
These restriction formulas were computed by Koike 
Thms. 11.2&11.3]. When A = A„, Koike's formu- 
las together with the degree filtrations yield 

M(2n + 1 ; A„) < 2;+ 1 = V (2n + 1 ; ski + i„) 



and 

M(2n;Xi] 



<2t+l 



■ 0y(2«;iA! + A„_ 1 )©V(2n;iA 1 +A„). 



6. Finally by the Weyl dimension formula, 

■? ~~r\ \ „ (2n — 1 + s 
dbmV{2n+\;sl\+Xn) =2" 



and 

dimV (2n;sAi +X„) = dimV (2n;sAi +A„) = 2 

Combining the dimension formulas yields 
dimM(n; A)<2r+i 



n-l 



2n — 2 + s 



F > 



as desired. 



dimV 



□ 

Remark. Corollarv l5.5l matches the Moller bound 1 18] for cu- 
bature on S", but it is more general because the measure 
need not be centrally symmetric. 
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Figure 7: A 175-point 30-cubature formula found by Wandzura and Xiao 1 31]. 



6. A LOCAL LOWER BOUND 

Our final application of moment maps is to help establish 
a local lower bound on the density of points of a PI or PB 
cubature formula on the simplex A„. The bound was origi- 
nally inspired by PI cubature formulas due to Wandzura and 
Xiao 1 3 1 ] which were found by simulated annealing. As in 
the example shown in Figure the points in these formulas 
accumulate transversely at the edges of the triangle. Another 
related result is that the limiting density of the points of Gauss- 
Legendre quadrature (i.e., the zeros of Legendre polynomials) 
is — J ||8]. This density can be interpreted as the linear 



projection of uniform measure on a circle, which is related to 
Archimedes' moment map (Figure[9}. 

Theorem l6.1l and Corollary 16. 2l establish a lower bound on 
the limiting density of any sequence of PI and PB formulas 
on A„ that generalizes the limiting density of Legendre zeros. 
Moreover, if the local density is high in certain regions, in par- 
ticular near the vertices of A„, then the weights there must be 
low. By this reasoning, Theorem l6.3l and Corollarv l6.4l estab- 
lish that a f-design on A„ requires many more points than an 



efficient f-cubature formula does) as t — > °° (namely 0(t 2 ") 
points versus 0(t") points). Along the way, Theorem 16.31 
establishes that Gaussian quadrature for an arbitrary weight 
function is very sharply locally optimal among all positive 
quadrature formulas. Finally Scholium l6.5l generalizes the re- 
sults for uniform measure on A„ to uniform measure on an 
arbitrary simple convex polytope. 

Theorem 6.1. A PI or PB 2t — 1-cubature formula F on the 
simplex A„ is an e-net with respect to the metric 



ds L = 



dx\ 







2xq 2xi 



dx\ 
2jc n 



in barycentric coordinates, where cos2e is the highest zero of 
the Jacobi polynomial p\ n 1 (x). 



In the proof and later, we will abbreviate (n 
superscripts. 



1,0) as "#"in 



Proof. The idea of the proof is to find, for each p 6 A„ and 
each e' > e, a P of degree 2t — 1 on A„ such that 



Ja„ 



P(x)dx > 0, 
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Figure 8: A polynomial on A n with a small positive island. 



but P(x) > only when x € A„ is in the e'-ball Z? £ / (/?) around 
p. We can call this ball the positive island of P(x); see Fig- 
ure [8] The existence of such a polynomial P forces F to have 
an evaluation point in B £ (p), for otherwise P(F) < 0. 

We first claim that the stated metric is the distance between 
fibers of the moment map % with respect to the Fubini-Study 
metric on CP". To see this it suffices to check the following: 
The real locus RP" C CP" is perpendicular to the fibers of % 
and meets each generic fiber 2" times. Indeed, n is a bijection 
on the orthant MP" with non-negative projective coordinates. 
Moreover, RP" is isometric to the orthant 5" of a unit n- 
sphere, and the restriction of % agrees with the restriction of 
the Xu map %\. The metric ds z on A n is exactly the push- 
forward of the standard metric on S" under Ti. See Figure|9] 
for an example. 




Figure 9: The moment map n restricts to the Xu map %\ . 



Consider the linear projection a : A n — * [— 1 , 1] given by 

a(x) = 2x -l. (1) 
The map a sends uniform measure on CP" to the measure 

ll{x)=n2 l - n {\-x) n -\ 

The t th orthogonal polynomial with respect to this measure ji 
is the Jacobi polynomial Pf = P t . Let pj be its highest 
zero. 

Define a polynomial Q 5 : CP" -> K by 

P f # (x) 2 (x-pf + g) 

Gs (z) = Gs (*) = (x _ p # )2 (2) 



where x — a(7c(z)) and 8 > 0. It has degree 2t — 1 as a polyno- 
mial in jc, as well as a polynomial on CP". Moreover, Qg van- 
ishes at the zeros of Pf, except at the highest zero, at which its 
value is positive. Therefore by Gaussian quadrature (!) with 
respect to the measure ji, 

[ Qs{z)dz= C Q s (x)dn>0. 

J CP" J -I 

At the same time, Q is non-positive outside of the region 

x > pf -8. 

This region corresponds to the ball of radius e' around (1:0: 
: ■ ■ ■ : 0), with 

2(cose') 2 - 1 = cos2e' = p* - 8. 

This can be confirmed by comparing with the orthant KP" 
from mentioned previously. Note that e' — > e as 8 — > 0. 

Given q S CP", define Qg j by rotating Qg by some isom- 
etry of CP" that takes (1 : : : ••• : 0) to q. Define Q T sp : 
A„ — > R, where p = 7c(q), by averaging Q over torus fibers: 

Then 

/ (% p (Z)<a= [ Q s ,p(z)dz>0, 
Ja„ ' Jcp" 

and ^ is non-positive outside of the ball of radius e' around 
p = 7C(q) in the induced metric on A„. Thus, Q§ p has the 
desired properties. □ 

Remark. A somewhat weaker version of Theorem 16 . 1 1 holds 
when F is positive and exterior, but with real evaluation 
points. Polynomials similar to Q T 5 can be constructed directly 
as products of factors that vanish on quadratic surfaces in 
R" D A„, with only one unsquared factor that vanishes on the 
boundary of B £ i(p). As it happens, the boundary of B E i(p) is 
a quadratic surface. We did not refine this sketched argument 
into a proof with explicit estimates. 

Corollary 6.2. Any sequence of PI or PB t-cubature formulas 

— 1/2 

on A n has limiting point density £l(t" Ylk x k )> where x € A„ 
is fixed and given in barycentric coordinates, and t — > °° 

Proof. The corollary follows from computing the volume 
form corresponding to the metric ds 2 in the statement of The- 
orem and estimating the covering radius e. The asymp- 
totic behavior of zeros of Jacobi polynomials is given in 
Abramowitz and Stegun 1 1, p. 787]. The key step in the esti- 
mate is the limit 

P ( "' fo) fcos^) 
lim ; =2 a e- a aU a (6), (3) 

where J a (z) is the ordinary Bessel function of the first kind. 
Convergence to the limit is analytic in 9. Thus 
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for every fixed k, where cos f ? /t is the &th zero of P^ a ' 6 ) (x) 
and is the kth zero of 7 a (x). The estimate can be estab- 
lished directly in our geometry by noting that P, # (2|zo| 2 — 1) is 
a harmonic function on CP" . The harmonic equation on CP" 
is then locally approximated by the harmonic (or Helmholtz) 
equation on R 2 ", whose radial solutions are derived from 
Bessel functions. 
In our case, 



cos2e = cos — -, 
t 



for (2t— 1) cubature. So 



e = ^(l + (l)) = e(r 1 ), 



which is also 0(f ) for f-cubature 



□ 



Theorem 6.3. Let /I be an arbitrary normalized measure on 
R whose support has at least 2t points. Let p\,..,,p t and 
W\, . . . ,Wt be the points and weights of Gaussian t-quadrature 
for the measure jX. Let F be a positive t-quadrature formula 
for ju. Then for each 1 < k < t, F has at least one point in 
the half-open interval (pk-i,Pk]> where po = — °°. Moreover, 
the total weight of all points in (— °°,/>i] is at most w\, with 
equality if and only ifF is the Gaussian quadrature formula. 

Note that Theorem 16. 3 1 is further sharpened by symmetry: 
F must also have at least one point in each half-open interval 
[pk,Pk+i), with p t+ \ = oo, and its total weight in [pt,°°) is at 
most Wt. 

Proof. Let <p t (x) be the t th orthonormal polynomial with re- 
spect to ji (with either sign), and let A, be the leading coeffi- 
cient of (j), (x). If k — 1, let 



P(x) = 



<Pt(x) 2 (pi+8i -x) 
(x-pi) 2 



with 5i >0. If k> 1, let 



P(x) 



<j>t (x) 2 (x -pk-i-5o) (p k + 8 1 -x) 



(x- p k -\) 2 {x- p k ) 2 
with 8\ ^> 5q > 0. In both cases, 

P{x)dn > 

by Gaussian quadrature, while P is only positive on the inter- 
val (pk-i + 8o,Pk + 8\). Therefore F has at least one point in 
this interval. Since F only has finitely many points, the limit 
8i — > establishes that F has a point in (pk-i, Pk] ■ 
For the second claim, let 



P(x) 



<M*) 2 

(x-p t 



\2- 



Then by Gaussian quadrature, 



P{x)dn=w x P{p l ). 



Let q\,...,qk be the points of F which are at most p\, and let 
Vi , . . . , Vfc be their weights. Then 

k k 

P(x)dfi = P(F) > £ Wj P( qj ) > P^) £ wj. 

j=i ;=i 

The first inequality holds because P is non-negative; the sec- 
ond because P decreases on ( — °o , p i ] . □ 

Corollary 6.4. The least weight of any positive t-cubature 
formula on A n (with uniform measure) is 0{t 2n ), uniformly 
in t. Any t-design on A„ has £l(t 2n ) points. 

Proof. If F is a f-cubature formula on A„, the map a (see 
equation Q) sends it to a f-quadrature formula a(F) on 
[—1,1] with Jacobi-polynomial measure. If F is positive, then 
the least weight of a(F) is at least that of F. On the other 
hand, Theorem l6.3l establishes that the least weight <x(F) is at 
least the last Christoffel weight w t . 

The first claim follows by estimating this weight. One of 
the standard formulas for the general Christoffel weight Wk is 

A t+1 \\<h(x)\\l 



w k ■ 



M!(pk)<Pt+i(pk)'' 



where <p t (x) is the f th orthogonal polynomial, A t is its leading 
coefficient, and pk is its kth root. In our case, ty, = Pf, pk = 
pj k , and k = t . We compute: 



\P?\\l 



2t + n 
_,(n-\+2t 
t 



0(2'). 



To estimate (ff)'(ptt) anc ^ u+l(Ptt)> we a g a in appeal to the 
limit in equation Q- Differentiating both sides by 9, we ob- 
tain 

(pfy (cos e/t) (sin e/t) 



lim 



tPf{\) 



al a Q- a a\j' a (d). 



(t+n-\\ 



Note that P*(l) 
the various parts of the limit yield 



= &(t n ). For a fixed value of 0, 



e 



(/f)'(cos 7 ) = 0(* 



n+U 



By the same token 



e, 



(Pf)'(cos-)=0(? B+1 ) 

when t approaches a fixed value of 9, as is the case when 
0, = 0f t is given by 

p tt =cos — . 
By a similar calculation, 

{pf + M) = ®{t n - 2 ). 

The conclusion is that 

w, = 0(r 2 "), 

as desired. □ 
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Scholium 6.5. Let K C W be a convex n-polytope with N 
facets. Let F be a t-cubature formula on K with uniform mea- 
sure. If F is PI or PB and if K is simple, then F is an e-net 
with respect to the metric 

X\ X2 x N 

where £ = (9(1/?). If F is positive, then its least weight is 
0(t~ 2n ). If it is a t-design, then it has at least 0(t 2n ) points. 



x 




Figure 10: A simplex L 1 (A„) whose facets contain facets of K that 
meet at x. 



Proof. (Sketch) The proof of Theorem l6. ll retains its strength 
if K C A„ and we pass from A„ to K, provided that the positive 
island of the polynomial Q T S - lies within K. In this case 

/ Qs,p(x)dx> f Q T s /x)dx>0. 

In order to properly position Q T 5 - for all p S /T, we need sev- 
eral embeddings of K into A„. For each vertex x G K, choose a 
linear embedding L that sends x to some vertex of A„, and that 
sends the facets incident to x to facets of A„. (Equivalently, 
for each vertex x € K, choose a simplex (A„) D K whose 
facets includes all facets of K that meet at x. See Figure [Toll 
Then there exists a finite set of L such that the positive islands 
of polynomials of the form Qs vq ° ^ together cover K. The 
formula F must have a point in each island, which establishes 
that F is an e-net. 

Similarly, the proofs of Theorem l6.3l and Corollary I6.4l re- 
tain their strength if uniform measure on K projects by a map 
a to a measure v on [—1,1] which is dominated by 

H(x) =2- n n{\-x) n - 1 

and which agrees with ji in a neighborhood of 1 . (Of course 
ji cannot dominate v if V is normalized, so this condition on 
V must be dropped.) In this case 

/ P(x)dv > [ P(x)dfi > 

Jr Jr 

for the first half of Theorem 16. 31 for the interval (—°°,pi), 
while 

/ P(x)dv < / P{x)djX 
Jr Jr 

for the second half of Theorem l6.3l A suitable projection a 
can be realized by positioning K in A„ so that it touches the 
vertex xq = 1, and then restricting the usual map a to K. □ 



7. OTHER COMMENTS 

In this article we have studied the toric moment map on 
CP", and on C" restricted to S 2 "- 1 (which can be interpreted 
as the level surface of an invariant Hamiltonian on C") as it 
applies to the cubature problem. Many of the constructions 
apply equally well to arbitrary toric varieties. To begin with, 
every complex projective variety X inherits both a metric and 
an affine real structure from CP". If X is toric, it also has a 
volume-preserving moment map whose image is a centrally 
symmetric polytope. However, the variety X rarely has much 
more symmetry than its moment map image. 

The duality between toric cubature (in particular trigono- 
metric cubature) and lattice packings explored in Section [3] 
suggests a different limit of the cubature problem. Let K C W 
be a centrally symmetric convex body. For simplicity let 
F = {p a } be a periodic discrete subset of M" with a periodic 
weight function p a i— > w a - Since it is periodic, it has a well- 
defined Fourier transform F. In this setting, F is a Fourier 
/T-cubature formula if and only if F is disjoint from the in- 
terior of K. The (continuous) Fourier cubature problem is to 
minimize the density of F among all /T-cubature formulas or 
all positive /T-cubature formulas. If F is a lattice with equal 
weights, then F has the same property and Fourier /T-cubature 
problem reduces to finding the best lattice packing of K. It 
would be interesting to find examples of non-lattice formulas 
that are better than the best lattice formula. 

We conjecture that a version of Corollary 15 .51 holds, using 
Theorem 15.21 and the same spinor bundles, for any centrally 
symmetric subvariety X C S". That is, we conjecture Moller's 
bound for these varieties, even when the measure jj. is not cen- 
trally symmetric. 

Theorem l6. ll shows why some tempting approaches to con- 
struct efficient PI or PB formulas on the simplex A„, even the 
triangle A2, are bound to fail. For example, if the points of 
a putative cubature formula F are fixed in advance, the ques- 
tion of whether it admits non-negative weights for f-cubature 
reduces to linear programming. But if the points are arranged 
in some lattice with spacing l/k, Theorem 16.11 shows that 
the weights can only be non-negative if k = H(f 2 ), so that 

\F\ = a(t 2 »). 

We believe that the requirement that K be simple in 
Scholium 16. 51 is not essential. More generally we conjecture 
that similar results hold if K is not convex. We also conjec- 
ture that the bounds in Theorem 16.11 and Scholium 16.51 are 
sharp to within a constant factor. The cubature formulas found 
by Wandzura and Xiao support this conjecture, at least when 
K = A„. 

The proof of Theorem l6.1l was partly inspired by the linear 
programming method to bound kissing numbers, f-designs, 
and sphere packings 0, Q3 0, H2]. Xu observed that 
the method for f-designs also yields bounds on PI f-cubature 
J34ll . In fact it yields an upper bound on the £2 norm of the 
weights of a PI f-cubature formula, which implies a lower 
bound on the number of points. We conjecture that linear pro- 
gramming methods could be used to improve the constants in 
Theorem| 6.1l 

Krylov 1 16] established that if {F t } is a sequence of interior 
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f -cubature formulas for a measure jj., then {f(F t )} converges 
to J x f(x)dji for every continuous / if and only if the l\ norm 
of the weights of F t is bounded as t — > «>. We conjecture then 
that Theorem l6. II still holds assuming a bound on the l\ norm 
of the coefficients of F instead of assuming that F is positive. 
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